Vorticity and Magnetic Field Generation from Initial Anisotropy in Ultrarelativistic 

Gamma-Ray Burst Blastwaves 
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Because conical segments of quasispherical ultrarelativistic blastwaves are causally disconnected 
on angular scales larger than the blastwave inverse Lorentz factor, astrophysical blastwaves can 
sustain initial anisotropy, imprinted by the process that drives the explosion, while they remain 
relativistic. We show that initial angular energy fluctuations in ultrarelativistic blastwaves imply 
a production of vorticity in the blastwave, and calculate the vortical energy production rate. In 
gamma-ray burst (GRB) afterglows, the number of vortical eddy turnovers as the shocked fluid 
crosses the blastwave shell is about unity for marginally nonlinear anisotropy. Thus the anisotropy 
must be nonlinear to explain the magnetic energy density inferred in measured GRB spectra. 
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Special relativistic shocks occur where the gravity of a 
compact object, such as a black hole, accelerates plasma 
to ultrarelativistic velocities. Emission form such shocks 
is observed in astrophysical sources on many length, 
time, and energy scales, and is often attributed to the 
synchrotron process. The composition of the plasma, 
the magnitude of the preshock magnetic field, and the 
Lorentz factor of the shock are generally unknown. In 
the best-studied system — the radiative afterglow that fol- 
lows the bright flash of 7-rays in gamma-ray burst (GRB) 
phenomenon — the structure of the emitting region can 
be modeled from the observed emission. In the standard 
GRB afterglow model [3, [3] , the radiation is produced in 
a relativistic blastwave shell propagating into a weakly 
magnetized plasma. The afterglow emission is then the 
synchrotron radiation from nonthermal electrons gyrat- 
ing in a strong magnetic field of the shocked plasma. 
Detailed studies of GRB spectra and light curves have 
shown that the magnetic energy density in the emitting 
region is a fraction €3 ~ 1CP 2 to 1CP 3 of the internal 
energy density Q and that the magnetic field must be 
present over at least a few percent of the blastwave thick- 
ness [3| . The origin of this downstream magnetic field is 
a longstanding important open question. 

Compressional amplification of the weak pre-existing 
microgauss magnetic field of the circumburst medium 
merely yields es ~ 1CP 9 0]. The strong small-scale field 
generated in collisionless plasma instabilities in the shock 
transition was considered as a candidate for the post- 
shock field Q. However, simulations of shock transition 
layers Q indicate that the small-scale field decays rapidly 
over few plasma skin depths (A s ), and does not persist on 
distances from the shock transition (~ 10 9 A S ) where the 
emission originates. Suggestions have been made that a 
persistent magnetic field develops in the shock precur- 
sor due to the streaming of shock-accelerated protons [8| 
or cascade-generated e ± pairs [j| . The feasibility of these 
scenarios depends on the insufficiently understood details 



of particle acceleration and streaming instabilities. 

Recently, Sironi h Goodman [ItJ showed that if the un- 
shocked medium is strongly inhomogeneous, significant 
vorticity is produced in the shock transition. The re- 
sulting downstream turbulence amplifies the weak seed 
magnetic field of the preshock medium to the observed 
level. The density of the interstellar medium may vary at 
the required level in GRB associated with the iron core 
collapse in mass-losing very massive stars (i.e., in "long"- 
type GRBs), but is not expected in the GRB associated 
with old progenitors (i.e., in "short" -type GRBs). While 
the strength of the magnetic field in the afterglows of 
short-type GRBs is not well-constrained, the similarity 
of the afterglow light curves and spectra in the long- and 
short-type GRBs [2| suggests that the magnetic field is 
amplified by the same process in both systems. 

Here we propose a source of vortical energy in the shock 
downstream that could generically be present in any ul- 
trarelativistic blastwave: the initial angular anisotropy of 
the energy carried by the blastwave. A unique feature of 
ultrarelativistic blastwaves is that they are composed of 
many small, causally-disconnected patches. The patches 
lose causal contact when the driving outflow (e.g., an 
electromagnetic jet launched near a black hole) accel- 
erates the swept-up ambient medium, and evolve inde- 
pendently until the blastwave decelerates sufficiently (see 
below). In GRBs, an additional source of angular inho- 
mogeneities is the phase in which the initial ( "prompt" ) 
7-rays are emitted, when a sizable fraction of the outflow 
energy is dissipated and emitted as 7-rays. This takes 
place after the relativistic outflow achieves its maximal 
Lorentz factor, and thus any variation in the dissipation 
in causally disconnected regions translates directly into 
angular energy fluctuations in the shock wave that ulti- 
mately blasts into the ambient medium. 

There is also direct observational motivation for initial 
angular anisotropy in GRB outflows. Significant angular 
fluctuations have been invoked to explain the large vari- 
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ation of the 7-ray luminosity between bursts [lfj and the 



intraburst variability observed in many afterglows 11|. 
The afterglow polarization (l2T | indicates a breaking of ax- 
ial symmetry; the correlatedness of this polarization with 
large amplitude variability suggests blastwave anisotropy 
as the source of variability 13|, [14J. E.g., the afterglow 
in GRB 021004 exhibits a strong light curve variability 
correlated with a variable linear polarization; anisotropy 
on scales ~ 2w/l with I ~ 200 and a nonlinear energy 
contrast ~ 3 can explain the data 14 1 . 

We proceed to estimate the conditions for successful 
turbulent magnetic field amplification in quasispherical 
ultrarelativistic blastwaves with initial angular energy 
fluctuations. Consider an ultrarelativistic point explo- 
sion with total isotropic-equivalent energy E propagat- 
ing into a medium of uniform density po- We work in 
the rest frame of the unshocked fluid and the explosion 
center. In spherical symmetry, the Lorentz factor of the 
strong shock wave that forms at the leading edge of the 
blastwave decays in time t as T = (j^E/^po) 1 / 2 and its 
position is located at R — t(l — |r ) (the speed of light 
is unity). The Lorentz factor in the shock downstream 
7 = r/(2x) 1//2 and the pressure p = |r 2 /x 17/ ' 12 are ex- 
pressed in terms of the Blandford-McKee [l5[ variable 
X(r,t) = (t- r)/(t -R) + 0(T- 2 ). Here, r is the radius 
from the center of the explosion and T 3> 1 is assumed. 

Consider initial (t — > 0) linear fluctuations in the en- 
ergy of the explosion, E{9,4>) = [1 + Se(9)]E, where at 
no loss of generality we assume axial symmetry, d<j, — > 
and 5e{0) = 5eYi (8). Here Y[ m (9,4>) is the spheri- 
cal harmonic; the case with an azimuthal dependence, 
m =/= 0, can be treated equivalently. Conical angular 
segments of the blastwave separated by A9 > §r -1 , 
which corresponds to I < 37IT, have not been in causal 
contact since the beginning of the explosion, and are 
evolving independently as fragments of spherical explo- 
sions with energy fluences displaced from the spherical 
average. Following [l6j| . we expand the radial Lorentz 
factor j(r,9,t) = j(r, t)[l + <5 7 (r, t)Yi (9)], the pressure 
p(r,9,t) = p(r 7 t)[l + S p (r,t)Yi (9)], the shock radius 
R(0,t) = R(t)[l + S R (t)Y l0 (9)], and the velocity of the 
fluid in the 9 direction u(r,9,t) = u(r,t)deYio(9), which 
is here assumed to be Newtonian, as linear perturba- 
tions around their spherical averages. Then, at times 
t » t crit = (^E/Ppo) 1 / 3 , i.e., when A9 > fL' 1 
and for x( r i^) ~ L t ne fluctuations simply inherit their 
initial values and are independent of time and radius, 

= ~T2$e, 5r = ^r _2 <5E, <5 7 = 0, and u = 0. 

Gruzinov (l6j linearizes the equations of relativistic 
hydrodynamics, d a T a ^ — 0, where T a @ is the energy- 
momentum tensor, and changes variables (r, t) — > (£,t), 
where £, = j hix(r, t) and r = — | mT(i), to derive equa- 
tions governing the evolution of the fluctuations 
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5 p + 35 p -8Si f + —S 7 -2l(l + l)u = 0, 

u + u'-^-u+^e 3T +^S p = 0, (1) 

where the overdot and the prime denote derivatives with 
respect to r and £, respectively, and the equations are 
valid for £ <C —jT an d t < 0, as the blastwave is New- 
tonian at larger £ or t. At the shock transition (£ = 0), 
the continuity of the energy-momentum flow across the 
transition implies the shock jump conditions [l6| 
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The first two of Eq. Q] can be expressed in terms of the 
Riemann variables f± = 2(5 7 ± propagating along 

the C ± characteristics with velocities = 3 =F 2y/S 



/* + c ± (/ ± ) / - + /") + C*« = 0, (3) 
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where a ± = \ T -j= and C ± = (1 T + 1). Following 

[l(| , we adopt the notation f = f + and g = /_ . 

Dynamics of the shocked fluid governed by Eq.Q]is eas- 
ily understood: an initial energy-pressure fluctuation 5 P 
sources transverse (5 P , u) and longitudinal (<5 P , <5 7 ) acous- 
tic fluctuations, which oscillate after the mode comes 
within the causal horizon at t = t cr i t (<5 7 equals the ra- 
dial velocity perturbation in the local fluid rest frame). 
This is analogous to the horizon re-entry of superhorizon 
cosmological fluctuations in the early universe. 

The initial value problem (Eq.[T]) subject to the bound- 
ary conditions (Eq. [2]) is well-posed if a boundary con- 
dition at the origin of the C + characteristics, at some 
distance £ max from the shock is specified. Gruzinov 
[l6| reports that the solution near the shock is indepen- 
dent of the boundary condition if £ max is large enough. 
The independence can be understood as follows. Since 
a + /c + > 0, the fluctuations / + are damped from large 
to small £; this is because the explosion energy is con- 
centrated near the shock [lij]. The reverse /~ fluctua- 
tions grow much slower than the / + decay, so the former 
do not seed the latter at large £. Having specified a 
boundary condition, e.g., /(£ m ax) = 0, the evolution of 
perturbations in the shock downstream is uniquely deter- 
mined and can be solved for, e.g., after Taylor expansion 
ip(€> T ) = E^LofaO^^nCr), where ip = f,g,u. 

An approximate solution is obtained after assuming 
that the perturbations decay after they start to oscil- 
late, and that the decay is complete everywhere ex- 
cept within few wavelengths A^A + from the shock, where 
A + = 2it\c + \/k + , and n + is the frequency of oscillations 
of fluctuations propagating along the C + characteristic. 
The frequency is estimated by substituting the last of 
Eq. [T] in the first of Eq. [U which yields k+ ps [±(1 - 

+ l)] 1 /2 e 3r/2^ Thus we get ^ NX+ K e -3r/2. 

the decay proceeds from large to small distances from the 
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shock. Expanding to the linear order yields 

fo + c + f 1 -a + (f +go) + C + u Q = 
go + c~gi - c\T (/o + go) + C~ u o = 
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where the first three equations are the original equa- 
tions of motions for /, g, and u, and the last three 
are the two constraints imposed at £ = (the shock 
jump conditions), and a constraint imposed at £ = £ max . 
Integrating the ODEs in Eq. 2] with initial conditions 
/o(— oo ) = — go(— °°) and uq(— oo) = and iV = 1, 
we find that the resulting approximate solution <5 P (£ = 
0, r) = -L[/ (t) — go( T )] closely matches the solution in 
(l6| . thus validating our approximations. 

We employ the definition of relativistic vorticity uj = 
V x H 17], which is conserved, dtUJ — V x (w x u5) = 0, in 
the smooth (i.e., shock- free) part of the flow in an ideal, 
barotropic fluid with ultrarelativistic equation of state. 
Here, H = p 1 / 4 7?/, where as before, p is the fluid pres- 
sure, 7 = (1 — v 2 — u 2 )~ x / 2 is its Lorentz factor, and we 
have decomposed fluid velocity into v = wf + Since 
ii « d, 7 k (1 — w 2 )" 1 / 2 is the radial Lorentz factor. Vor- 
ticity vanishes in an unperturbed, spherically-symmetric 
blastwave, but in a blastwave with initial anisotropy, it 
is produced at the shock and advected into the down- 
stream. In the (£,t) coordinates, the vorticity conserva- 
tion equation becomes Co + uj' + 9uj + 0(L~ 2 ) = 0, where 
uj is the amplitude of the sole nonvanishing component of 
the vorticity, uj — udeYio4>. The vorticity decreases into 



the downstream, k>(£, r) = e 9 ^uj(0,t — £). In terms of 



the fluid perturbation the vorticity equals 
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where we have kept only the leading-order terms. The 
vorticity in Eq. [5] is conserved if the fluid perturbations 
evolve according to Eq. Q] Because 5 p (t — > — oo) ^ 
while (5 7 ~ and u ~ 0, vorticity does not vanish identi- 
cally prior to the horizon entry of the mode. 

We follow the prescription in [ItJ to calculate the vor- 
tical energy, which is defined in the local, noninertial rest 
frame of the shocked fluid; we denote the quantities in 
this frame with tilde, f = j(r — vt) and t = j(t — vr). 
Vorticity transforms as a 2-form f2 = uj^dr Ad(r9). Since 
dr A d{r6) —> 7 _1 df A d(r6) to the leading order in 7 _1 
(terms involving dt do not contribute to vorticity), we 
have uj ~ w/7, and the ^-dependence of the vorticity is 
given by w(f , t) = e -( n / 2 )«(i(0, r - £). 
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FIG. 1: Evolution of the pressure perturbation 8 P , radial 
Lorentz factor perturbation <5 7 , transverse velocity u, and the 
fluid- frame vorticity uj, and the fractional energy in vortical 
motions at the shock. The quantities 8 P , 5 7 , and uj are ex- 
pressed in the units of the initial pressure anisotropy 8 Pl o, and 
e VO rt is expressed in the units of 5 p0 . 



The fluid is subject to differential acceleration, and 
thus there is not a single inertial frame in which we can 
calculate the vortical energy. As an approximation, we 
work in the instantaneous inertial rest frame of the fluid 
element at £ = 0. Simultaneity in this frame is equivalent 
to di(£, t) = 0, which implies that dr w ^T^tdC This 
is not exact; the numerical coefficient depends on the 
approximation, and thus the forthcoming estimates are 
crude. We also replace covariant derivatives with normal 
derivatives in what follows. We carry out approximate 
projection of the vorticity onto the shock plane 
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«K = 0). (6) 



The vortical energy density in the local fluid rest frame 

equals [17] e vort = ^p\H so \\ 2 p~ 1/2 , where H so i is the 

solenoidal component of H, and p is the proper energy 
density. The solenoidal component can be written as 

a curl of a vector potential (H) so \ = V x A. Setting 

V • A = 0, the vorticity is related to the potential via 
V 2 A = —uj. Recall that Q oc d$Yio, which is not an 
eigenfunction of the angular component of the Laplacian. 
However, for I ^S> 1, at no loss of generality we can restrict 
analysis to a small patch of the spherical shell, e.g., near 
the equator, \0 — \ \ < Z — 1 , where Yiq(0) w — cos(W) and 
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the shell is locally quasi-planar. Then A — A sin(Z0)</> 
and A = ^tae~ lr / 1 is the solution of (we freely set the 
origin of the coordinate f = at £ = 0) 

^4- l ^A=L Q{ rJ)^LaS(P), (7) 

Or z t 7T 7T 

where S(f) is the 5-function. Close to the shock we take 



curl to obtain ((H SI 



(^) 2 ~ (sk^h) 2 ^ where 



the angular brackets denote an average over cos 9. Substi- 
tuting in e vort and diving by the internal energy p yields 
the fractional vortical energy e vort = £ V ort / P 
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where K = | (g^) 2 ~ 2.5 x 10~ 4 is a numerical constant 
sensitive to the specifics of our approximations. 

Fig. [TJ shows the evolution of the fluid perturbation 
variables, the fluid- frame vorticity and the fractional en- 
ergy in vortical motions. The latter oscillates and reaches 
peak average amplitude e VO rt ~ 3 x 10~ 2 (5 Pi o 2 after an l- 
mode becomes causal, for T ~ (0.01 — 0.1) x I. Here, 
6 P fl = S P (T = oo ) is the initial fractional pressure fluctu- 
ation. The plots in Fig. [T] are scaled to be independent 
of I, but abscissal range satisfying the premise that the 
blastwave is ultrarelativistic does depend on I. 

A fluid with nonvanishing vorticity develops eddies 
which give rise to magnetic field amplification via the 
turbulent dynamo mechanism [hsl ]. Strong amplifica- 
tion is expected when the number of eddy turnovers is 
larger than unity. The number is estimated as the ra- 
tio iVcddy ~ tcrossAeddy) of the time scale t cmss ~ \t/l 
on which the fluid crosses the shell of the blastwave to 



the eddy turnover time i c ddy 
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±1(1 + 1), we have 
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For e V ort ~ 0.03, we expect iV 0( jdy ~ 1- Therefore, within 
the limits of our approximations, a nonlinear initial per- 
turbation, Spfi > 1, is necessary to amplify the magnetic 
field. If amplification happens, the final energy density 
in the field will be in approximate equipartition with the 
vortical motions, es ~ e V ort- With 5 Pi q ~ 1 we find that 
tB ~ 3 x 10 -2 , consistent with the typical values inferred 
from GRB spectra. The anisotropy of GRB blastwaves 
could be fully nonlinear. Since our treatment is inade- 
quate in the nonlinear regime (see [19i | for an analytic 
approach to aspherical blastwaves), numerical simula- 
tions of relativistic blastwave hydrodynamics are needed 
to study the nonlinear turbulence produced in the shock. 

Fig. [1] shows that the vorticity associated with an l- 
mode persists through an increase by xlO 3 in T 2 , which 
corresponds to an increase by x 10 4 in the observer time 



(t bs oc r~ 8 / 3 ). Thus a single mode I ~ 500 produces vor- 
ticity during the observed GRB afterglow (10 2 — 10 6 s af- 
ter the burst). Since r max ~ 100, we indeed expect such 
an I mode from causality considerations. Our estimates 
indicate that Se ~ 1 is required for strong magnetic field 
amplification, although in this regime, our linear treat- 
ment is not strictly applicable. A measurable signature 
of such angular anisotropy is afterglow flux variability 
at the level ~ Se on time scales i bs for which I ~ 27rr 
[llj ]. At later times the variability amplitude decays as 
~ 2tt{T/1)Se and can drop below 1% while the vortical 
energy is still high (l/T ~ 100). The observed variability 
depends on the initial spectrum of the fluctuations. In 
many GRB afterglows a variability is observed as pre- 
dicted in our model. In some cases no flux variability is 
observed around 10 5 s after the burst down to the level 
1% — 10%. There, as we explain above, large-Z fluctu- 
ations may be responsible for magnetic field generation 
without imprinting an observable variability. 
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